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Abstract 

We analyze the quantization of the Pais-Uhlenbeck fourth order oscillator 
within the framework of deformation quantization. Our approach exploits 
the Noether symmetries of the system by proposing integrals of motion as the 
variables to obtain a solution to the *-genvalue equation, namely the Wigner 
function. We also obtain, by means of a quantum canonical transformation 
the wave function associated to the Schrodinger equation of the system. We 
show that unitary evolution of the system is guaranteed by means of the 
quantum canonical transformation and via the properties of the constructed 
Wigner function, even in the so called equal frequency limit of the model, in 
agreement with recent results. 
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1. Introduction 

Whenever one consider curvature terms, for example in general relativity 
or brane inspired models, one is faced naturally with held theories described 
by Lagrangians with higher order derivative terms. The Pais-Uhlenbeck 
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fourth order linear oscillator, originally introduced in [l| , is perhaps the sim¬ 
plest example and dehnitely the best known higher derivative mechanical 
system and, in particular, it has served as a toy model to understand several 
important issues related to Ostrc^radsky instabilities emerging naturally in 
higher order held theories [21, y, 14, la, la, IT], Isj . Recently, the Pais-Uhlenbeck os¬ 
cillator has been used as a guide to study higher order structures associated to 
supersymmetric held theory [^, PT-symmetric Hamiltonian mechanics 10 


and geometric models within the scalar held cosmology context 0 . In this 
sense, it is important to mention that naive quantization procedures for the 
Pais-Uhlenbeck model has to be enhanced in order to recover unitarity in 
a physical allowed sector. Our main motivation is thus to address for the 
Pais-Uhlenbeck oscillator, within the perspective of the deformation quanti¬ 
zation formalism, the long standing problems associated to the non-unitarity 
of higher derivative theories. As we will see, within our formulation the uni¬ 
tarity is guaranteed straightforwardly, even in the equal frequency limit of 
the model, by the introduction of a well-dehned Wigner distribution. 

The framework of deformation quantization was introduced in [l^ as an 
alternative approach to the problem of quantization. In this formalism one 
uses, as guidelines, the Dirac quantization rules in order to pass from classical 
physics to the quantum realm. As is to be expected, a consistency require¬ 
ment for such a quantum theory is the existence of a classical limit, that 
is, a quantum system should reduce to its classical counterpart whenever 
the limit of h, the Planck constant, tends to zero. From this perspective, 
the quantization of a classical system could be seen as a deformation of the 
algebraic structures involved in a parameter encoding the quantum nature 
associated to the system {h in our case). Furthermore, the quantization rules 
require that for any classical observable there is a corresponding quantum 
observable, and similarly, that the Poisson bracket corresponds to the quan¬ 
tum commutator. All these requirements can be achieved by replacing the 
usual product of the algebra of smooth functions on the classical phase space 
with an associative non-commutative product, depending on h, such that the 
resulting commutator is a deformation of the Poisson bracket. We refer the 
reader to [l^, where results on the explicit construction of maps between 
classical and quantum observables are explained in detail, to Refs. [l4, 15|, 
where conditions on the existence of the star product are exposed, and to 
the reviews [l^, [13, H for general aspects of deformation quantization, as 
well as for more recent developments. 

Our approach is based on taking advantage of the symmetries inherent to 
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the Pais-Uhlenbeck model in order to construct the Wigner function that 
contains the relevant quantum information of the system. In this man¬ 
ner, we show that there exists a couple of integrals of motion associated 
to Noether charges, which in turn serve as privileged variables in order to 
hnd the solutions to the T*r-genvalue equation. Further, in order to obtain the 
quantum wave functions we consider both, classical and quantum canonical 
transformations. At a classical level we transform, in a standard way, the 
Pais-Uhlenbeck system to a simpler model composed of the difference of two 
uncoupled harmonic oscillators for which the Wigner function may be also 
obtained. We then use the latter Wigner function to obtain a wave equation 
which, by means of a quantum canonical transformation, may be used in 
order to obtain the quantum wave function for the original Pais-Uhlenbeck 
system. The resulting wave function is identical to the one obtained in by 
a different reasoning. We also show that in the equal frequency limit of the 
model the source of the continuous spectrum can be traced out through a lin¬ 
ear canonical transformation that maps the Pais-Uhlenbeck Hamiltonian to 
a Hamiltonian composed of a discrete spectrum part plus a continuous spec¬ 
trum part, contrary to the unequal frequency case. Besides, we demonstrate 
that in the equal frequency limit the Wigner function is certainly unitary as 
consequence of composition of unitary transformations considered through 
the quantum canonical transformations. In this sense, our results explicitly 
manifest the ghost-free feature of the Pais-Uhlenbeck model, in complete 
agreement with j3, |^. 

The article is organized as follows. In Section [2l we include a brief review 
of deformation quantization in order to set our notation and to define some 
useful structures. We also consider quantum canonical transformations as 
they will be essential in our context to obtain the wave functions associated to 
the Pais-Uhlenbeck oscillator. In Section [31 we analyze the Wigner function 
for our model in terms of its integrals of motion and we identify the quantum 
wave equation. Also, in this section we detail the equal frequency limit for the 
Pais-Uhlenbeck oscillator. In Section IH we include some concluding remarks. 
Finally, we address some technical issues related to the construction of the 
Pais-Uhlenbeck wave function in the Appendix. 
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2. Deformation quantization 

2.1. Basic notions 


in 


The origins of deformation qnantization were first introdnced by H. Weyl 
19| . The main idea is to associate a qnantnm operator W[f] in the Hilbert 
space to every classical observable f{q,p) dehned on the phase space 

This operator is known as the Weyl operator and it is explicitly given 
by 


f{q,p) ^ W[f] := 


1 

{2TihY^ 


(Tp d^^exp 


h 


{q-p + p-i) 




( 1 ) 


where / is the Fourier transform of / G given by 


hv,o 


(Tq d>exp 



{p-q + i-p) 


f{q,p), 


( 2 ) 


q, p are operators satisfying the canonical commutations relations, and the 


integral is understood in the weak operator topology [20, l2l|]. Later on, 
Wigner obtained an inverse formula which maps a quantum operator 
into its symbol, that is, a differential operator with polynomial coefficients 
dehned on classical phase space. This map, known as the Wigner function, 
results a quasi-probability distribution function, explicitly written as 


= '^"y r ( 9 - 5 !/) e-‘yp,p (q + . ( 3 ) 

Among the most important properties of the Wigner function we would like 
to point out the following. In the hrst place, Wigner function is a particular 
representation of the density matrix which is normalized, bounded and real. 
In the second place, for a given quantum wave function 'i/j{x), Wigner function 
represents a generating function for all spatial auto-correlations. Finally, 
in the classical limit h i-G- 0 Wigner function reduces to a highly localized 
probability density in the coordinate space j^. 

Subsequently, Moyal found an explicit formula for the symbol of a quan¬ 
tum commutator between operators [^ . 


if *9){q,p) = / exp 


2 ^ q'- ^ Pi 



( 4 ) 
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This ^-product is known as the Moyal product, and will be dehned as an 
involution in ([7]). For a two dimensional case (n = 1) this deformed product 
reads 

1 ‘ ^ \ ^ ^ Jc\ 

^ ' m=0 ^ ' 

Straightforwardly, one may check that the Moyal product satisfies the fol¬ 
lowing properties: 

1. f g = fg + 0{h), that is, the T*r-product is a deformation of the usual 
pointwise product between smooth functions. 

2. With respect to the T*r-product, the Weyl operator hF is a homomor¬ 
phism of algebras, 

W : ^ (£(L2(M^)),o) , (6) 

where 

W[f^g]:=W[f]oW[g], (7) 


between distributional functions defined on the phase space and the 
space of bounded operators on the square integrable functions 
under the composition o as a product. 

3. The deformed Poisson bracket associated to the T*r-product is given by 

[/, gh-=^{f^g-g^f) = {/, g} + 0{h ), ( 8 ) 


where {•, •} denotes the standard Poisson bracket, and [•, is known as 
the Moyal commutator or the Moyal bracket. As a consequence of (??), 
Moyal bracket may be interpreted as a Lie bracket on the space of Weyl 
operators (isl ]. 


In order to avoid convergence problems, the Moyal product is defined 
not in the space of smooth functions, but in the extended space, 

C'°°(R^")[[h]], which corresponds to the space of formal power series in h with 
coefficients in or alternatively, it can be defined on the space of 

Schwartz functions using an inteyal representation and then extend it to a 
suitable space of distributions j^, The interpretation of this product as 
a non-commutative deformation on the algebra of observables was introduced 
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by means of Gerstenhaber’s algebras. In particular, the existence of 


a T*r-product for an arbitrary symplectic manifold was demonstrated in [15 


Subsequently, as a consequence of the formality theorem, Kontsevich solved 
the existence and classification problem of star products on a generic finite 


dimensional Poisson manifold 28 


2.2. Time evolution and canonical transformations 

The expectation values of observables and the time evolution of states 
can be computed in a similar manner as in classical mechanics, however, 
the usual pointwise product between functions and the Poisson bracket are 
replaced with the ^-product and the Moyal commutator, respectively. In 
consequence, the expectation value of an observable A in a state ip, is given 
by the expression 




♦ (27rft)»/2 


r‘qd'‘pA{q,p)*p(q,p), 


( 9 ) 


where p{x,p) is the Wigner function dehned in ([3]). In this sense, expectation 
values of physical observables in deformation quantization correspond to ob¬ 
taining the trace of an operator with the density matrix, in close analogy to 
standard probability theory [^ . 

The dynamical equation of the quantum distribution p{x,p), results in 
the counterpart of Liouville’s theorem in classical mechanics describing the 
time evolution of a classical distribution function, and is given by the formula 


dp{q,p) 

dt 




H -k p — p-k H 
ih 


( 10 ) 


where is a distinguished real function from the algebra of observables, 
namely, the Hamiltonian. Being time-dependent, this evolution equation, 
also known as Moyal’s equation, does not completely determine the Wigner 
function for a system |29j. Then, just as in the conventional formulation of 
quantum mechanics, a systematic solution may be inferred from the spectrum 
of the stationary problem. Static Wigner functions obey a more suggestive 


*-genvalue equation, inducing Bopp shifts [16 


H{q,p)*p{q,p) = p{q,p)-^H{q,p) = H (^q+^-^dp,p-^-^dg^ p{q,p) 

= Ep{q,p), (11) 
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where E corresponds to the energy eigenvalue associated to the Hamiltonian, 
leading to the spectral properties of the Wigner function as a quantum distri¬ 
bution function jl^. These quantum properties are related to the fact that 
Wigner function is not positive semi-definite, allowing in principle negative 
values for certain areas of the phase space. This counter intuitive negative- 
probability aspect has been speculated as a way to detect quantum interfer¬ 
ence which may be measured and reconstructed indirectly in the laboratory 


As within the framework of deformation quantization the observables are 
represented by smooth functions dehned on phase space, and thus they trans¬ 
form classically, the outcome of a canonical transformation on the quantum 
*-genvalue equations results in an appropriate transformation of the Wigner 
function (3ll |. 

A general transformation of the phase space coordinates is defined as a 
smooth bijective map T : 3 {q,p) {Q,P) £ which transforms 

every observable A e by 


A' = A oT. 


( 12 ) 


Further, the star product in the new variables V satisfies the natural condition 

3 

{f^g)oT={foT)^'{goT), f,ge , (13) 

in such a way that the T*r'-product is given by 

/ V ^ = / exp - y g , (14) 

where the vector helds Dqi, Dp. correspond to the transformed derivations 
dqi, (9p. according to the rule stated in equation flT^ . 

Among the phase space coordinate maps there are some which play a 
special role in classical mechanics, namely the canonical transformations or 
symplectomorphisms, that is, those maps which preserve the symplectic form, 
and thus the Poisson bracket structure 


{'/‘.pq = {Q‘./’q='55. (15) 

In phase space quantum mechanics, a quantum canonical transformation is 
a transformation T such that it preserves the form of the corresponding 
deformed Poisson bracket, namely, the Moyal’s bracket 

= ( 16 ) 
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where [•, denotes the Moyal’s bracket transformed by T to the new coor¬ 
dinate system [32|, l33| , 


lf,gl. = [foT-\goT-^]^oT. 


(17) 


The transformation of coordinates T induces a unitary operator on the 
Hilbert space Ut '■ —)■ transforming vector states and ob¬ 

servables to the new coordinate system in such a way that for any quantum 
observable A the transformation reads 


A'{Q,P) = UTA(q,p)U^\ 


(18) 


This expression allows us to derive the form of the operator JJt, which cor¬ 


responds to the celebrated Dirac’s quantum transformation [34 


{UT^l^){Q) = 


(27rh)"-/2 


f / 

d^F , 

1 d^Qj 

dqdQ 




(19) 


being F{q, Q) the generator of the classical canonical transformation such 
that 


P = 


d F{q,Q) 
dq 


P = 


dnq.Q) 

dQ 


( 20 ) 


In this expression, the implementation of the quantum canonical transforma¬ 
tion on state eigenfunctions is realized through a generalization of a Fourier- 
type transformation containing ^-corrections with respect to the classical 
canonical counterpart, thus encompassing the quantum nature of the opera¬ 
tor. In addition, it can be demonstrated that Ut is indeed a unitary operator 
on the Hilbert space L‘^(EA) j^, 321. 


3. Deformation quantization of the Pais—Uhlenbeck Oscillator 

In this section, we perform a deformation quantization of the Pais-Uhlenbeck 
oscillator by solving the corresponding i*r-genvalue equations and, in order to 
obtain the associated wave functions, we make use of the theory of quantum 
canonical transformations discussed in the previous section. 

















3.1. Wigner function 

In order to study the quantum dynamics of the system, we will start by 
considering the Lagrangian 

L{q, q, q) = + ^l)q^ + , ( 21 ) 

where g is a real-valued function of time, q and q its hrst and second time- 
derivatives, respectively, and the parameters hli and VL 2 correspond to a pair 
of frequencies which are taken real and positive. More explicitly, this model 
is characterized by the following fourth-order differential equation of motion 


dS 


+ + Ql)^ + QlQlq 


0 , 


( 22 ) 


which can be derived directly from the Lagrangian. The canonical Hamil¬ 
tonian can be obtained by means of the Ostrogradsky’s method 


35 . To 


implement this method, the phase space involves, in addition to the canon¬ 
ical coordinates {q,Pg), an extra canonical pair of variables, namely x := q 
with corresponding canonical momentum 0,0, a. The momenta are de- 
hned a.s px := dL/dq and Pq := dL/dq — d/dt (dL/dq), and in our case they 
explicitly read 


Px = q, 

Pq = -{nl + nl)q-^, (23) 

respectively, while the canonical Hamiltonian is obtained through the Leg¬ 
endre transformation H{q,pq, x,Px) = Pqd + Pxi — L{q, q, g), and thus 


H{q,pq,x,px) 


Px 

PqX+ — + 




(24) 


For this theory, it is natural to dehne the generalized Poisson bracket 

dA dB dA dB dA dB dA dB 

+ 


{A 5} := 


dq dpq dpq dq 
in order to obtain canonical relations 


dx dpx dpx dx ’ 


(25) 


{9.P,} = 1 = , 


(26) 
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while the rest of the brackets among the phase space variables are vanishing. 
This Poisson structure enable us to write the canonical Hamilton equations 
of motion which, by making the identification x = q, straightforwardly lead 
to equation fl2^ . Therefore, the Hamiltonian H, obtained in describes 
the dynamics in phase space for the fourth order Pais-Uhlenbeck model, and 
thus this Hamiltonian corresponds to the one which is to be quantized, using 
in our case, the techniques related to deformation quantization as reviewed 
above. 

Our goal now is to calculate the Wigner function. As we will see below, 
the *-genvalue problem (ITT]) can be solved directly without first solving the 
corresponding Schrodinger equation. Considering the Pais-Uhlenbeck oscil¬ 
lator with Hamiltonian given by 0241) . the T*r-genvalue equation OTT]) for the 
Wigner function p{q,Pq,x,px) explicitly reads 




where we have used the associative ^-product 



(28) 


•k := exp 


which is naturally associated to the canonical Poisson structure fl2^ of the 
Pais-Uhlenbeck oscillator. The eigenvalue equation may be decomposed 
into two partial differential equations corresponding to the real and imaginary 
parts of the equation, respectively. These equations explicitly read 



(29) 


for the real part, while 


xOq Pxdx T T ^ 2 ) Pi.Q y Pgy ^y Px:') 0 ) (^0) 
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corresponds to the imaginary part. Although, these equations may seem 
challenging at hrst sight, we can take advantage of the deeply relation be¬ 
tween Moyal quantization and the canonical structure of the theory in order 
to hnd solutions to these equations. By using the symmetries of the Hamil¬ 
tonian, and employing the conserved Noether charges as variables for the 
Wigner function, one can infer the solution to the T*r-genvalue problem (j^ . 
Indeed, using Noether theorem, suitable generalized to higher derivative the¬ 
ories, one hnds that the Pais-Uhlenbeck oscillator remains invariant under 
the symmetry 


37 


q + e 


d?q 

dfi 


± - ^^ 2 ) 


dq 

dt 


(31) 


which, as a consequence, imply the existence of two global integrals of motion 


Ji = 


J') — 


1 


{px + + ipq + 




ni-ni 

1 

ni-ni L 


(pq -\- ^2 {px + 


2W2 


(32) 

(33) 


These integrals of motion are associated to inhnitesimal transformations of 
the form fIM]) . However, when one considers inhnitesimal transformations 
along the velocity vector held, g, together with inhnitesimal time transfor¬ 
mations, the integrals of motion are associated to the so-called energies for a 
second order diherential Lagrangian. In particular, one may check that the 
Pais-Uhlenbeck Hamiltonian fl2T)) is related to the integrals of motion of our 
interest by the identity H = (Ji — J2)/2. A general statement of Noether 
theorem and the energies for a second order Lagrangian system is reviewed 


in [38|, Appendix A], while a detailed exposition may be found in [39 . 

Thus, using constants of motion fl3^ and ([33)), the Wigner function given 
by 


Pnmiqy Pqy X, Px) ^ 


m+n 


TT" 


n? 


iw) (S) ■ 


corresponds to a solution to the i*r-genvalue equation ([27|). Here, the L’s 
stand for the Laguerre polynomials dehned by the Rodrigues formula 


Ln{z) = 


ezd^(e-^z^) 


n\ 


(35) 
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As we will see in the next section, Wigner function (IM|) is intrinsically re¬ 
lated to the Wigner function of the harmonic oscillator. Besides, it turns out 
that the Hamiltonian of the Pais-Uhlenbeck oscillator is canonically equiv¬ 
alent to the Hamiltonian corresponding to the difference of two uncoupled 
harmonic oscillators. This imply, as long as the canonical transformation 
remains linear 33|, that the Wigner function of the Pais-Uhlenbeck is equiv¬ 
alent to the Wigner function of a pair of harmonic oscillators, but written in 
terms of its own Noether charges. 

Also, after some algebraic manipulations, one may check by substituting 
the Wigner formula (l34|) into equation ([29]), that the system has energy -k- 
genvalues 


E-nm ^ ^ 2 / T 2 ) ^2 ) 


n, m = 0,1, 2, 


(36) 


as it is well-known for the Pais-Uhlenbeck oscillator. 


3.2. Wave function 

In this section we are interested in obtaining the wave function associated 
to the quantum mechanical system of the Pais-Uhlenbeck oscillator within 
the framework of deformation quantization. To this end, we will use the 
formalism of quantum canonical transformations already discussed in sec¬ 
tion 12.21 This can be achieved by following the next strategy. First, by 
means of an appropriate canonical transformation, we will map the Hamilto¬ 
nian of the Pais-Uhlenbeck oscillator into a Hamiltonian describing a pair of 
uncoupled harmonic oscillators. The second step will be to solve the Wigner 
function corresponding to the new Hamiltonian, and then, through a Fourier 
transformation, we will calculate the wave function related to the harmonic 
oscillator problem. Finally, we will consider a quantum canonical transforma¬ 
tion in order to obtain the wave function belonging to the original problem, 
that is, the Pais-Uhlenbeck oscillator. The aim of this approach is twofold. 
On the one hand, instead of calculating the wave function directly from the 
Pais-Uhlenbeck Wigner distribution, we note that the quantum canonical 
transformation not only leads to more manageable integrals, but also it al¬ 
lows us to compare with previous results found in the literature where the 
wave function is obtained by solving the Schrodinger equation. On the other 
hand, the equal frequency limit Oi —)■ 02 is analyzed via canonical transfor¬ 
mations at both, the classical and the quantum levels j^, 1^. 
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Consider for definiteness, the case of nneqnal frequencies and assume 
> r 22 . The Hamiltonian fl24p can be brought into diagonal form by 
applying the canonical transformation j^, 


Q 


Px 


1 - Pi 

^^1 ’ 

fiiPi - nlX2 

^ynj - ’ 


- P 2 

’ 


Pq 


^ - nl 


which is realized by the generating function 


F(g, X, Xi,X 2 ) = fli'yqX 2 + jxXi - Vt\qx - VtiXiX 2 , 


(37) 


(38) 


where 7 := — Using such a canonical transformation, one hnds 

that the Hamiltonian is mapped into the difference of an uncoupled pair of 
harmonic oscillators 


H(Xi,Pi,X2,P2) 


Pi + nlxl Pi + nlxl 

2 2 


(39) 


Following the T*r-genvalue equation ffTTl) for this Hamiltonian, the resulting 
equation is 




-{P 2 - 


■^2 ( ^2 + ~^^ P 2 


-2E 


p(Xi,X2,Pi,P2) = 0. 


By virtue of its imaginary part 

h {—Pidxi + VL\Xidp^ + P2dx2 ~ ^\X2dp^') P = 0, 


(40) 


(41) 


the Wigner function p is seen to depend on two variables, Zi = APl^^^/Mli and 
Z 2 = 4.H1^^/hVL2, where FTf" ;= (P^ + Vl\XI)/2 and — (7f + ^i ^|)/2 
correspond to the Noether charges to the Hamiltonian of uncoupled harmonic 
oscillators ([39]). Then, the *-genvalue equation, reduces to a couple of simple 
ordinary differential equations 


— Z'3‘^ — 3 — 


E 

2hVti 


P{zi,Z2) = 0 , 


i = 1,2 


(42) 
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This equation acquires the same form as that for the simple harmonic oscil¬ 


lator [Ig, therefore, one may be easily convinced that the Wigner function 


is determined by 


-Pi;-^2;-P2) — 


(-1) 


n-\-m 2 H? 


hCl 




1 e 


4^os( 

Mil 


4^osc 

h^lo 


.(43) 


where the energy spectrum E results equal to the energy spectrum of the 
Pais-Uhlenbeck oscillator fl36|) . 

We can obtain the wave function 'ijj{Xi,X 2 ) in the coordinates {Xi,X 2 ) 
by Fourier transforming in the momentum variables the Wigner function fl43p 
adapted to our system, obtaining jiol ] 


^(^1,^2) 


1 

-^*(0, 0) 


dPidP2 p 



jPiXi+P2X2/ri 

^ 5 


(44) 


where the constant ■^*(0,0) may be determined up to a phase by normal¬ 
ization of '0(^1; ^2)- By using the following identity between Laguerre and 
Hermite polynomials, denoted by FT’s (4lj |. 



Hn{x 


a)H4x + 


2'^y/nn\e (2(a^-|-6^)) ,(45) 


the integral gives the wave function in the position space 


'4^nm (^1 ; ^2 ) 


_ 1^72 H-m 


q 2 j 5^2 


(46) 


The final step consists in calculate the wave function associated to the 
Pais-Uhlenbeck oscillator. To this end, we make use of the quantum canoni¬ 
cal transformation operator, defined in ffT9|) . where the generator, F(g, x, Xi, X 2 ), 
of the classical canonical transformation takes the specific form of ([38]). Ex¬ 
plicitly, the wave function for the Pais-Uhlenbeck oscillator in the position 
space is given by the formula 


i^nmiQyX) = N JdXidX 2 |exp 

t/^nm (Wi, W2) } , 


- {ni^qX2 + -fxXi - njqx - ^1X1X2) 


(47) 


where iV is a normalization constant. From now on, we will consider h = 1 for 
simplicity. Also, here we only state our results, leaving the technical details 
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on the construction of the wave equation to the Appendix. Bearing this in 
mind, after performing the integration in fH7|) . we are able to determine the 
wave function for the Pais-Uhlenbeck oscillator, which explicitly reads 


A 




x) = Nnm exp [-iVtiVt 2 qx\ exp 

where A := fii — ^ 2 ; and the functions (pnm stand for 

m\{n — m)\ 


(pnm{q,x), (48) 




k=0 


(m — k)\k\{n — m + k)\ 




+ 


I n—m+k'*''*' k 


Hf^ , m <n, 


4^nm(^qy x'j ^ ^ A 


n\{m — n)\ 


k=0 


n — k)\k\{m — n + k)\ 


, m>n. (49) 


Here, the constant A a ■= iA/di/HTn^ only depends on the frequencies, 
while the H~^^s and the represent Hermite polynomial evaluated in the 
arguments 


H+ := 


H- := 


i^/^l (^ 2 ? — ix) 
{klgQ + ix) 


(50) 


respectively, and the constants Nnm behave as normalization factors. The 
wave function fHHll stand for the solution to the Schrodinger equation asso¬ 
ciated to the Hamiltonian operator of the Pais-Uhlenbeck oscillator [^. As 
it was proven in j^, these solutions are normalized resulting in a pure point 
spectrum, and all eigenfunctions form a complete orthogonal basis in the 
Hilbert space L^(M^). 


3.3. Equal frequency limit 

In the equal frequency limit, A = fli — f22 0, the quantum wave 

function (I48|) cease to be normalizable, implying that the spectrum acquires 
continuous values. Indeed, within our context the Wigner function (l34)) oscil¬ 
lates wildly but eventually approximating to zero, therefore, the probability 
distribution appears as a generalized function, in such a way that the energy 
of the system i s eq ually likely to be found anywhere in an interval [E, E + dE] 


for any E 1^. Explicitly, whenever A = 0 the source of the continuous 
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spectrum can be traced out through a quantum canonical transformation 
m defined via the linear classical canonical transformation 

qQ 2 ^qx QxQi Q 1 Q 2 


F{q,x,Qi,Q2) = 


Vtqx VtxQi 
\/2 4 ^ ^2 


where we have considered the single frequency as hi := hli = 172 
generating function maps maps the Pais-Uhlenbeck Hamiltonian 
new Hamiltonian 


Q 2 

Hn = n{Q^P2-Q2Pi)-^{Q\ 


Ql 


(51) 

This 
to a 


(52) 


which clearly differs from that obtained as the difference of two uncoupled 
harmonic oscillators in (|39ll . We also may easily deduce that the spectrum 
of the equal frequency Hamiltonian (15^ is composed of a discrete spectrum 
coming from the angular momentum part Q 1 P 2 — Q 2 P 1 , and of a continuous 
spectrum originated from the squared norm of the position variables Ql + Ql 
as demonstrated in Ref. j^, 


P^k = ^h(m- . ( 53 ) 

Moreover, as it was indicated in 1^, the Pais-Uhlenbeck system results 
ghost-free, as even in the equal frequency limit A —)■ 0 the evolution operator 
is certainly unitary. Again, within our context, the unitary property is readily 
obtained from the Wigner functions fIMD and fITH]) . as both are related by 
a quantum canonical transformation of the form flTSD . where the unitary 
operator Ut is in fact given by the integral transformation fl47D . Following 
this reasoning, since the time evolution of the pair of harmonic oscillators is 
a solution of Moyal’s equation [l^ 


p{Xi,X2, Pi, P 2 , t) 


U-^^piXi,X2,Pi,P2,0)^U., 


(54) 


where the ^-evolution operator 

U^{Xi,X2,Pi,P2) = h^h+--- 

= (55) 

n,m 

corresponds to a unitary operator. Therefore, the evolution of the Pais- 
Uhlenbeck oscillator is also unitary, as a consequence. 
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4. Discussion 

In this article, we analyzed the fourth order Pais-Uhlenbeck model within 
the Deformation quantization formalism. We obtained both, the Wigner dis¬ 
tribution function and the wave equation, for this system. In particular, in 
order to obtain the Wigner function we explicitly consider the symmetries 
of the system associated to its Noether charges which were established as ad 
hoc variables for solving the T*r-genvalue equation. For the wave equation, we 
proceed by first considering a classical canonical transformation that map the 
Pais-Uhlenbeck model to a system of uncoupled harmonic oscillators. After¬ 
ward, our strategy was to obtain the Wigner function for the new system, 
and then to construct the associated wave function for this system. Finally, 
we reached the wave function for the Pais-Uhlenbeck model by a quantum 
canonical transformation. This wave equation resulted identical to the one 
considered in . We also showed that the model contains a continuous spec¬ 
trum and results ghost-free, the two conditions together being a consequence 
of the unitariness of the relevant quantum canonical transformations. 

Further studies are necessary in order to see the generality of conserved 
quantities associated to Noether charges as appropriate variables to find so¬ 
lutions to the *-genvalue equation for a generic model. This will be done 
elsewhere. 
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Appendix A. Technical steps for the construction of the Pais- 
Uhlenbeck wave function 


In the following we will detail some calculations that were omitted in 
the main text. Explicitly, by means of the quantum canonical transforma¬ 
tion flT7|) . the wave function for the Pais-Uhlenbeck system fH7|) can be ex¬ 
pressed as 




exp [—iillqxj 


j (iXi(iX2exp 


i^xX2 


^2X1 

2 


iVti'yqXi 


17 




-iVtiXiX2 


VtiXl 


(^/n^X^) Hm (^^ 2 X 2 ) . (A.l) 


Defining Zi := ^/TtlXi and Z 2 := \/n 2 -A 2 , this expression may be written 


as 


V'(g,x) = 


exp [—iVilqx] 

\/DiD 2 


^ Z 1 Z 2 — 

D 2 2 


dZidZ 2 exp 




H^{Z^)H^{Z2) . 


(A.2) 


Multiplying both sides by using the generating func¬ 

tion for the Hermite polynomials 




exp -h 2tx) = ^ ’ 


fc =0 


(A.3) 


we perform the Zi and Z 2 integration using the Gaussian integral |41l 


/ dxexp [—p^x^ + qx) = exp 

-00 

thus obtaining in this way the identity 

00 00 , l , 

%L^ 

EE —= 27rexp [—iDiD25'3:] exp 


q VTT 

4p2 p 


A 


(A.4) 


k=0 1=0 


— -|- ^1^25'^) 


X 


X 


CXJ CXJ CXJ ^ 

EEEl^-i'Awj 

u—pt 7—n o—n V 

Hk{x^)Hi{x 


k=0 1=0 j=0 


if (Di-D2)^^+'^A \ 


kUljl 


^ ' ' (Dl + Q2fk+l+2j)/2 J 

(A.5) 


Using the double summation identities 

00 00 I 


cxD n 


k^n / ^ ^ ^k,n—k • 

n=0 k=0 n=0 k=0 

00 00 00 n 

E E - E E > 


(A.6) 


A:=0 n=k 


72=0 k = 0 


and relabeling the summation indices, we hnally obtain that the n, m-term 
of the sum flA.Sp . dehned as V'nm, results identical to the expression 
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